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Methods  are  presented  for  calculating  the  minimum  attainable  probability  of  error  in  detecting  a coherent 
quantum  signal  received  amid  thermal  noise  in  the  limits  of  low  and  high  signsl-to-noise  ratio.  In  the  Utter  limit, 
quantum-mechanical  perturbation  theory  is  applied  to  solving  the  detection  operator  equation  approximately. 
Graphic  results  are  furnished  for  binary  signals  transmitted  with  equal  prior  probabilities.  Previous  work  on  this 
problem  is  thus  extended  to  a broader  range  of  signal  and  noise  strengths.  The  results  are  applied  to  the  detection 
of  antipodal  two- photon-cohe rent-state  signals. 


UNCLASSIFIED 


SUMMARY 


OBJECTIVE 

The  objective  of  this  study  was  to  investigate  the  effect  of  thermal  background  noise 
on  optimum  quantum  receivers.  This  study  employed  the  methods  of  statistical  hypothesis 
testing,  estimation,  and  quantum  decision  theory. 


RESULTS 

Methods  are  presented  for  calculating  the  minimum  attainable  probability  of  error 
in  detecting  a coherent  quantum  signal,  received  amid  thermal  noise  in  the  limits  of  low  and 
high  signal-to-noise  ratio.  In  the  latter  limit,  quantum-mechanical  perturbation  theory  is 
applied  to  solving  the  detection  operator  equation  approximately.  Graphic  results  are 
furnished  for  binary  signals  transmitted  with  equal  prior  probabilities.  Previous  work  on 
this  problem  is  thus  extended  to  a broader  range  of  signal  and  noise  strengths.  The  results 
are  applied  to  the  detection  of  antipodal  two-photon-coherent-state  signals. 


RECOMMENDATIONS 

1.  Investigate  alternate  ways  of  solving  this  binary  communications  problem,  in 
order  to  bridge  the  gap  between  high  and  low  signal-to-noise  ratio  under  medium-to-high 
thermal  noise  conditions. 

2.  Investigate  the  effect,  on  thermal  noise,  of  the  detection  of  more  general  two- 
photon-coherent-state  signals. 
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INTRODUCTION 


The  rapidly  developing  area  of  electro-optical  (E/O)  technology  offers  great  potential 
for  application  to  military  command,  control,  and  communications  (C’>  systems.  Recent 
studies  have  shown  that  the  replacement  of  electronic  devices  with  E/O  components,  in  the 
information  transfer  process,  could  result  in  a significant  improvement  in  system  size,  weight, 
band  width,  covertness,  and  cost  It  is  likely  that  these  and  other  advantages  will  lead  to 
increased  use  of  electro-optical  technology  in  Navy  systems  of  the  future.  At  present,  a 
number  of  specific  naval  ship  and  aircraft  applications  (e.g.,  OPSATCOM  and  RED  PERCH) 
have  been  identified  and  active  research  and  development  efforts,  toward  implementing  this 
technology  in  these  roles,  have  been  initiated  at  NOSC  and  elsewhere. 

The  classical  detection  and  estimation  theories,  which  are  generally  utilized  in  optical 
system  analysis,  assume  that  information  carrying  fields  and  random  noise  can  be  measured 
with  arbitrary  accuracy.  Unfortunately,  the  propagation  and  detection  of  optical  signals  is 
actually  a very  complex  stochastic  process;  and,  it  must  be  treated  by  quantum  mechanics 
rather  than  classical  mechanics  to  be  rigorously  correct.  This  latter  fact  places  certain 
constraints  on  the  optical  measurement  process.  These  constraints  show  up  in  the  failure 
of  commutativity,  between  the  electric  and  magnetic  field  components,  at  different  space- 
time  points  on  the  same  light-cone.  Thus,  in  order  to  apply  the  principles  of  statistical 
hypothesis  testing  and  estimation  to  optical  signals,  the  theory  must  be  compatible  with  the 
laws  of  quantum  mechanics. 

When  we  speak  of  a quantum  “system”  in  connection  with  the  detection  and  esti- 
mation of  optical  signals,  we  refer  to  the  electromagnetic  field  at  the  aperture  A of  the 
observing  optical  instrument  during  the  interval  of  observation  (O.T).  This  field  can  be  de- 
composed into  spatio-temporal  modes,  ortho-normal  over  A and  (O.T);  and  for  these  modes, 
creation  and  annihilation  operators  obeying  the  usual  commutation  rules  can  be  defined. 
Thus,  the  establishment  of  the  Hilbert  space  of  state  vectors,  appropriate  for  the  application 
of  quantum  mechanics,  is  permitted.  In  many  problems  only  one,  or  a finite  number,  of 
these  modes  needs  to  be  taken  into  account. 

Unfortunately,  few  problems  in  quantum  detection  and  estimation  have  been 
solved  comprehensively.  This  paucity  of  solved  problems  has  limited  our  understanding  of 
the  import  of  the  new  theory  and  has  inhibited  the  development  of  intuitions  for  attacking 
more  complicated  problems.  In  detection,  for  instance,  nearly  all  solutions  apply  to  choices 
among  pure  quantum-mechanical  states  corresponding  to  the  absence  of  thermal  background 
light  In  most  of  these,  the  optimum  receiver  attains  error  probabilities  smaller,  by  an  order 
of  magnitude,  than  those  attained  by  conventional  detectors  (examples  are  given  in  ref  1 , 
ch  VI).  The  question  arises,  how  deleterious  is  a small  amount  of  thermal  noise?  This  is 
important  because  many  actual  external  noise  sources,  such  as  the  sun,  can  be  modelled  as 
thermal  blackbodies. 

1.  Helstrom  CW,  Quantum  Detection  and  Estimation  Theory,  chapt  V,  Academic  Press,  1976. 
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The  simplest  problem  in  ordinary  detection  theory  is  to  determine  the  optimum 
detector  of  a coherent  signal  of  known  form,  received  in  the  presence  of  Gaussian  noise, 
and  to  calculate  the  minimum  attainable  probability  of  error.  The  corresponding  problem 
in  quantum  detection  theory  has  not  been  solved  exactly.  Approximation  methods  have 
been  described  by  Yoshitani  (ref  2)  and  by  Klemm  (ref  3),  with  comments  by  Helstrom 
(ref  4);  and  the  problem  has  been  reviewed  in  a recent  book  (ref  1).  Here,  we  shall  extend 
previous  work  over  a broader  range  of  signal  and  noise  strengths.  The  principal  novelty  is 
the  application  of  quantum-mechanical  perturbation  theory  to  a suitably  reformed  detection 
operator  equation,  which  yields  the  asymptotic  form  of  the  error  probability  in  the  limit 
of  large  signal-to-noise  ratio  (SNR).  The  opposite  limit  of  low  SNR  is  also  treated.  Our 
results  permit  calculating  the  minimum  error  probabilities  attainable  in  detecting  antipodal 
two-photon-coherent-state  (TCS)  signals,  received  in  the  absence  of  background  noise. 

Before  proceeding  however,  a review  of  the  need  for  a quantum  mechanical  approach 
is  in  order. 


BACKGROUND 

In  classical  physics,  a system  can  be  characterized  by  performing  a series  of  dynamic 
variable  measurements  designed  to  completely  specify  the  state  of  this  sytem.  Each  measure- 
ment has  a particular  probability  density  function  associated  with  it.  Thus,  choices  among 
various  hypotheses  about  the  system  can  be  made  by  means  of  the  optimum  decision  stra- 
tegies prescribed  by  ordinary  statistical  hypothesis-testing  theory.  There  are  several  mono- 
graphs in  the  literature  which  describe  this  type  of  decision  process  in  some  detail. 

Nature,  on  the  other  hand,  precludes  such  a state  characterization  when  dealing  with 
a quantum  system.  In  particular,  the  propagation  and  detection  of  optical  signals  is  actually 
a very  complex  stochastic  process,  and  must  be  treated  by  quantum  mechanics  to  be  rigor- 
ously correct.  This  latter  fact,  places  certain  constraints  on  the  optical  measurement  process. 
These  constraints  show  up  in  the  failure  of  the  commutativity  between  the  electric  and  mag- 
netic field  components,  at  different  space-time  points  on  the  same  light-cone.  Hence,  in 
order  to  apply  the  principles  of  statistical  hypothesis-testing  to  optical  signals,  the  theory 
must  be  compatible  with  the  laws  of  quantum  mechanics.  Thus,  one  is  left  not  only  with 
the  problem  of  determining  the  best  way  to  process  the  results  of  one’s  measurements,  but 
also  with  the  problem  of  which  parameters  should  be  measured  in  the  first  place. 

As  an  example,  imagine  an  optical  communication  system  transmitting  messages 
that  have  been  coded  into  an  alphabet  of  M symbols.  The  message  source  emits  a new  symbol 
every  T seconds.  The  transmitter  can  produce  M different  optical  signals,  one  for  each  sym- 
bol in  the  alphabet.  They  are  modulations  of  the  electromagnetic  field,  produced  perhaps 
by  a laser,  and  last  no  longer  than  T seconds.  When  the  message  source  emits  the  jth  symbol 
of  the  alphabet,  the  transmitter  sends  out  the  j**1  signal. 

2.  Yoshitani  R,  On  the  Detectability  Limit  of  Coherent  Optical  Signals  in  Thermal  Radiation,  J Statis  Phys, 
vol  2,  no  4,  p 347-378,  1970. 

3.  Klemm  RA,  Detection  of  a Coherent  Signal  of  Known  Phase,  Info  and  Control,  vol  28,  p 35-S3,  May  1975. 

4.  Helstrom  CW,  Comments  on  Detection  of  a Coherent  Signal  of  Known  Phase  by  R.  A.  Klemm,  Info  and 
Control,  vol  35,  p 1 14-1 18,  October  1977. 
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The  receiver  collects  all  the  light  incident  on  its  aperture  A.  We  suppose  it  to  be  so 
synchronized  with  the  transmitter  that,  taking  account  of  the  time  needed  for  a light  signal 
to  pass  from  transmitter  to  receiver,  it  can  identify  the  beginning  and  end  of  each  signaling 
interval  of  duration  T seconds.  During  each  interval,  it  observes  the  incident  light  and  on 
the  basis  of  the  observation  decides  which  of  the  M signals  are  transmitted. 

A typical  observation  begins  with  filtering  the  light  to  cut  out  as  much  background 
radiation  as  possible.  Light  from  a local  laser  may  be  added  to  the  incident  beam,  as  to  a 
heterodyning  component,  and  the  combined  beam  then  focused  on  a photosensitive  surface. 
The  electrical  signal  so  generated  is  amplified  and  applied  to  a bank  of  M filters,  each 
appropriately  matched  to  one  of  the  M possible  signals.  Whichever  filter  has  the  largest 
output,  when  sampled  at  a certain  moment,  indicates  a decision  that  the  corresponding  signal 
was  the  one  transmitted  (ref  5).  Any  procedure  of  this  kind  is  subject  to  error,  because  of 
noise  created  by  the  background  and  by  the  physical  processes  of  the  receiving  system  itself. 
One  would  like  to  design  the  entire  procedure  so  that  the  average  probability  of  error  is 
minimal;  or,  when  costs  are  assigned  to  the  various  kinds  of  error,  the  average  cost  per 
decision  is  as  small  as  possible. 

Whatever  procedure  is  adopted  is  basically  a processing  of  the  electromagnetic  field 
at  the  aperture  A during  the  observation  interval  (O.T).  It  must  be  compatible  with  the  laws 
of  quantum  mechanics.  It  is  convenient,  as  suggested  by  Takahasi  (ref  6),  to  imagine  the 
incident  light  to  have  been  admitted  to  a lossless  box  or  cavity  behind  the  aperture  A,  which 
is  opened  at  t=0  and  closed  at  t=T.  The  field  within  the  cavity  is  identified  as  a quantum 
system  S.  Its  physical  properties,  which  determine  the  structure  of  its  Hilbert  Space  Hg, 
has  recently  been  described  (ref  1 ).  There,  one  finds  the  nature  of  the  signals  and  the  back- 
ground light  determining  a set  of  M density  operators  (pj,  p2, ....  p^j),  which  respectively 
specify  the  state  of  the  field  in  each  of  the  M possible  transmitted  signals.  It  is  necessary 
to  choose  one  of  these  density  operators  that,  in  some  sense,  best  represents  the  cavity  field. 
For  this  purpose,  the  field  is  observed  at  a time  T,  or  later,  by  an  instrument  I that  applies 
to  it  a probability-operator  measure  (p.o.m.)  (ref  7).  This  measure  will  have  M components, 
denoted  by  II  j , II2, ....  IIM. 

The  lossless  cavity  is  called  an  ideal  cavity.  The  artifice  of  observing  its  electromag- 
netic field  at  a single  instant  t,  T is  adopted  in  order  to  avoid  the  complication  of  analyzing 
quantum-mechanical  measurements  made  at  a succession  of  times  in  (O.T);  each  of  these 
would  interfere  with  the  field  in  an  unpredictable,  random  manner.  No  generality  is  lost 
for  the  cavity  field  at  t > T and  the  evolution  of  the  aperture  field  during  the  interval  (O.T) 
are  uniquely  related,  the  one  containing  all  the  information  provided  by  the  other.  More- 
over, no  scheme  for  processing  the  aperture  field  can  yield  a lower  average  cost  than  that 
attained  by  the  optimum  p.o.m.  applied  to  the  field  in  the  ideal  receiver.  In  order  to  show 
how  the  optimum  p.o.m.  is  determined,  we  consider  the  general  problem  of  deciding  among 
M density  operators  of  an  arbitrary  quantum  system. 

5.  Pratt  WK,  Laser  Communication  Systems,  Wiley,  1968. 

6.  Takahasi  H,  “Information  Theory  of  Quantum-Mechanical  Channels,”  Advances  in  Communication 
Systems,  ed  A Balakrishnan,  vol  1,  p 227-310,  Academic  Press,  1965. 

7.  Heistrom  CW,  Quantum  Detection  and  Estimation  Theory,  chapt  III,  Academic  Press,  1976. 
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There  are  M hypotheses  about  the  state  of  the  system  S,  of  which  the  jth  is  the  pro- 
position that  its  density  operator  is  py  In  our  example,  it  implies  that  the  jth  signal  was 
transmitted.  The  p.o.m.  affecting  the  decision  needs  only  M components,  II  j,  II2, ....  n^j, 
which  we  shall  call  detection  operators.  These  nonnegative-definite  Hermitian  operators  sum 
to  the  identity 

M 

S "i  - i.  <>) 

J=1 

and  specify,  according  to  reference  7,  the  conditional  probabilities 

Pr(j  |k)  = TKpknj),  Oik)  = 1,  2, ....  M,  (2) 

that  the  instrument  I chooses  hypothesis  Hj  when  is  true.  In  order  to  reach  that  decision, 
the  instrument  may  have  registered  certain  meter  or  counter  readings  and  processed  them  in 
some  way,  perhaps  by  a computer  that  ultimately  prints  out  one  of  the  integers  1,2,..., 

M.  These  detection  operators  Ilj  are  not  necessary  projectors,  but  rather  form  a generalized 
resolution  of  the  identity  (ref  7). 

An  instrument,  for  example,  that  simply  guesses  which  hypothesis  might  be  true, 
selecting  an  arbitrary  one  with  probability  I/M  without  even  interacting  with  the  system  S, 
in  effect  applies  the  p.o.m.  composed  of  the  operators 

nj=M_I  l,  j = I,  2, ....  M. 

These  are  certainly  not  projectors. 

Let  f;  denote  the  prior  probability  of  hypothesis  Hj  in  our  communication  system. 

This  is  the  relative  frequency  with  which  the  jth  signal  is  transmitted.  The  cost  of  choosing 
hypothesis  Hj,  when  Hj  is  true,  is  again  Cy.  Then,  the  average  cost  of  an  observational  strategy 
specified  by  the  p.o.m.  |na|  is 

_ M M M 

C=H  fjCijTr(pjni)=Tr  l Wjllj,  (3) 

i=l  j=l  i=l 

where  the  Hermitian  risk  operators  Wj  are  defined  by 

M 

wi=  I fjCijPj-  (4) 

j-1 

We  must  minimize  equation  (3)  under  the  constraints  that  equation  ( 1 ) hold  and  that  the 
detection  operators  Ilj  be  nonnegative  definite,  IIj>  0,  and  Hermitian. 

This  problem  resembles  one  of  linear  programming,  except  that  operators  are  involved 
instead  of  functions.  Its  structure  is  much  like  the  problem  of  minimizing  the  Bayes  cost  of 
a decision  strategy  in  ordinary  statistical  theory.  The  equations  for  the  solution  are, 

(Wi-T)ni*ni(Wi-T)  = 0,  i=  1,  2, . . . , M,  (5) 

Wj-T^O,  i = 1,  2, . . . , M,  (6) 


,C  - t *»  >•  * i v • . 


(5) 

(6) 


- 


with  the  operator 

M M 

T=I  njWj=I  Wjllj  (7) 

j=l  j=l 

required  to  be  Hermitian.  In  equation  (6)  the  inequality  sign  means  that  the  operator 
Wj  - T must  be  nonnegative  definite.  This  solution  to  the  problem  of  optimum  quantum 
hypothesis  testing  was  first  published  by  Holevo  (ref  8).  Equations  (5)  - (7)  were  reported 
by  Yuen  as  solving  a problem  in  vector-space  optimization  theory  later  recognized  as  the 
equivalent  (ref  9).  A demonstration,  based  on  the  principle  of  duality  in  vector-space 
optimization  theory,  can  be  found  in  a paper  by  Yuen,  Kennedy,  and  Lax  (ref  10). 


The  operator  T in  equations  (5)  — (6)  plays  the  role  of  a Lagrange  multiplier  taking 
account  of  the  constraint  of  equation  (1);  which,  involving  Hermitian  operators,  requires  T 
to  be  Hermitian.  We  call  T the  Lagrange  operator.  The  minimum  Bayes  cost  C_in  is 
(refs  1,  10): 

^min  = Tr  T . (8) 


BINARY  DETECTION  USING  THE  BAYES  CRITERION 

When  there  are  only  two  hypotheses  between  which  to  choose,  it  is  customary  to 
label  them  Hg  and  Hj.  The  system  S has  density  operators  pg  and  p j under  Hg  and  Hj, 
respectively,  and  their  prior  probabilities  are  fg  and  ? j,  ?g  + £ j = 1.  We  shall  see  that  the 
optimum  p.o.m.  for  deciding  between  the  two  hypotheses  is  projection-valued  and  can  be 
defined  in  terms  of  the  eigenvectors  of  a linear  combination  of  pg  and  p j. 

There  are  now  only  two  detection  operators,  rig  and  II  j , and 

n0  + ni  = l.  (9) 

From  this  condition, 

ngdig  + rij)  = (ng  + nj)rig 

yields 

n0ni=n,no;  (io) 

or,  the  detection  operators  commute.  The  Lagrange  operator  T is  now 

T = w0n0  + w1n1  (ID 

8.  Holevo  AS,  Statistical  Problems  in  Quantum  Physics,  Proc  2nd  Japan-USSR  Symp  Probability  Theory, 
Kyoto  vol  1 , p 2240,  August  1972. 

9.  MIT  Res  Lab  Electron  TR  482.  Communication  Theory  of  Quantum  Systems,  HP  Yuen,  30  August  1971. 

10.  Yuen,  HP,  Kennedy  RS,  and  Lax  M,  Optimum  Testing  of  Multiple  Hypotheses  in  Quantum  Detection 
Theory,  IEEE  Trans,  IT-21,  P 125-134 >1  arch  1975. 
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and 


f 

I 


I 

I 


w0-T  = w0-w0n0-w1n1=(w0-Wi)n1. 

Now,  equation  (5)  yields 

(W0-W1)n1n0  = 0; 

where,  by  equation  (4), 

W0-Wl  =fi  (Coi  -cii)pi  _?o(cio_coo)po 
= ^1  (C0i  ~ c 1 1 KP l 

with 

r0(C10-C00> 

fl(C0l -c,,)  • 


(12) 

(13) 


(14) 


(15) 


This  corresponds  to  the  classical  decision  level  Aq.  We  can  confine  our  operations  to  the 
subspace  of  the  Hilbert  space  spanned  by  the  eigenvectors  of  the  density  operators  pg 
and  p i,  corresponding  to  nonzero  eigenvalues.  That  subspace  may  coinc  ’ i with  the  entire 
space  3Cg.  In  it,  the  operator  Wq  - W j will  not  in  general  vanish,  and  therefore,  II  j rig  = 0; 
thus,  the  operators  Ilg  and  II  j are  projectors  onto  orthogonal  subspaces  of  3fg. 

The  nonnegativity  of  Wg  - T yields 

Wg  - T = (Wg  - w,)n,  = i ,(Cg,  - c, ,)(p,  - Xpg)n, > o . 


Since  Cgj  - Cj  j represents  the  net  cost  of  an  error  of  the  second  kind,  it  is  positive,  and 

n 1 (p  j - Xpg)  = (p  j — Xpg)II  j ^ 0 . (16) 

The  eigenvectors  | rjj ) and  the  eigenvalues  r\ j of  the  Hermitian  operator  pj  — Xpg 
are  defined  by 

(Pi  - Xpg)  1 TJj  > = J?j  t TJi  >.  (17) 

In  general,  some  of  the  eigenvalues  rjj  will  be  positive,  or  zero,  and  the  rest  negative.  For  each 
negative  eigenvalue,  form  by  equation  ( 1 6), 

<r?i  I (P i — Xpo>nl  It?i>  = i7i<r?i  lllj  lr?i  )>  0. 


With  tjj  < 0,  this  implies 
< I II  j I jjj  > <0, 


but  because  II  j must  be  nonnegative  definite,  this  requires 


IHj 

1 TJJ  > = 0, 

T?i<0. 

(18) 

Similarly, 

<*7i 

in0 

1 >7,  > = 0, 

»?i>0. 

(19) 
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Thus,  we  find  that  II  j projects  onto  the  subspace  of  Jf  s spanned  by  the  eigenvectors  I rjj  > 
corresponding  to  positive  eigenvalues;  and  Ilg  projects  onto  the  subspace  spanned  by  I tjj  > 
for  7?j  < 0.  The  subspace  spanned  by  any  eigenvectors  with  = 0 contributes  nothing  to 
the  average  cost.  We  can  write  the  optimum  detection  operators  as 

n0=S  U(-r?j)  IrjjXijj  I,  rij  = £ Uftjj)  IrjjXrjj  I,  (20) 

i i 

where,  U(x)  is  the  unit  step-function, 

1,  x>0, 

U(x)=  Vi,  x = 0,  (21) 

0,  x < 0. 

The  optimum  observation  strategy  can  be  described  as  the  “measurement,”  in  the  conven- 
tional sense,  of  the  operator  pj  — Xpq  (ref  11).  If  a positive  eigenvalue  results,  choose  H j ; 
if  a negative  one,  choose  Hq.  If  there  exists  a zero  eigenvalue,  and  it  turns  up  in  the  measure- 
ment, either  hypothesis  may  be  selected.  We  have  arbitrarily  picked  a strategy  that  chooses 
randomly  between  them  and  with  equal  probabilities.  The  minimum  Bayes  cost  is 

Cmin  = TrT  = foC00  + flC01  _^1^C01  - C1  rijUfrjj)  . (22) 

i 

Unfortunately,  finding  the  eigenvalues  and  eigenvectors  of  the  operator  pj  - Xpg  may  be 
very  difficult. 


BINARY  DECISIONS  IN  ZERO  THERMAL  NOISE 

Given  two  pure  states  Ii//q>  and  we  have  that 
Pq=  I^oX^qI 
and 

pj  = l^jXvi/j  I. 

As  this  is  a two  dimensional  space,  we  also  have 
Ilj  = tajXcoj  I 

with  j = 0, 1.  The  closure  or  completeness  relation  tells  us  that 
1 

1=  2 iwjX^  i = n0  + n, . 

j=0 

If  we  now  expand  k = 0,  1,  in  terms  of  the  lull’s,  we  obtain 

11.  Helstrom  CW,  Detection  Theory  and  Quantum  Mechanics,  Info  and  Control,  vol  10,  p 254-291, 
March  1967. 
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l'/'0>  = Xo0lw0>  + Xl0lwl  > 


(27) 


and 


1^1  > = Xoil"0>  + xllla,l>>  (28) 

where  x«  = <Wj  l^j  > . Let  us  associate  Iooq)  with  the  negative  eigenvalue  t?q  of  (p  j - Xpg) 
and  koj  ) with  the  positive  eigenvalue  rjj.  From 

(pj  - Xpy) b?  > = rj\rt  >,  we  have 

1\1/ Itj  > — XI\|/qXi|/o^  ) = ) • (29) 

Taking  the  scalar  product  of  equation  (29)  with  <i//qI  and  then  with  tyj,  yields 

W> -(r?  + X)<  =0 


(30) 


and 


(1-t?)  < ^Itj ) - X7  < \pQ  b?  > = 0 (31) 

with  7 = ( ^ j Ii^q  > . The  above  simultaneous  equations  will  have  a non-zero  solution  if 

= 0 


-0 J + X) 


| (1-T?)  -7^ 

This  implies  that 


Vi  _(1-X)+x/(1-X)2  + 4X(1-1yI2) \ 0 


and 


_ _ ( 1— X)  - >/(l-X)2  + 4X(l-lyl2)  <.  n 

r?0  2 

Therefore,  we  can  write 


Thus, 


and 


with 


(l-i7k)  -X7 


< t j I »?k> = Ak  7* 


<^0lT?k)  = (,-T?k)  Ak 


lAkr2  = 2R(r?k-g) 


■<  ^ 1 1 J?k  >' 

0 

<^0,T?k> 

0 

(32) 


(33) 


(34) 


(35) 


(36) 
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R - |[^(l-X)]2  + Xg)^ 
g = I-I7I2 
for  k * 0, 1 . 

The  false-alarm  and  detection  probabilities  are  equal  to 
Q0  = Pr(H1lH0[=Tr(p,n0) 

_ (l-r?j)2  [i?i  +X*g] 

X2I y 1 2 R 

and 

Qd  = Pr)H1lH,}  =Tr(p,n,) 

[»?j  + X*g] 

” R 

respectively. 


r 


COHERENT  DETECTION  IN  THE  PRESENCE  OF  THERMAL  NOISE 


In  the  preceding  subsection  we  discussed  the  decision  between  two  pure  states  in 
the  absence  of  any  noise.  Unfortunately,  this  is  not  very  realistic  since  most  optical  systems 
must  contend  with  some  sort  of  noise  mechanism  affecting  the  decision  process.  The  intent 
of  these  sections  is  to  investigate  binary  “on-off’  coherent  and  two-photon  coherent  detec- 
tion in  the  presence  of  thermal  noise.  In  particular,  we  will  expand  the  approach  developed 
above  to  take  into  account  the  effects  of  a thermal  noise  source  in  the  decision  process. 


DETECTION  OF  BINARY  ON-OFF  SIGNALS 

In  an  on-off  binary  communication  system,  a coherent  optical  signal,  such  as  the 
output  of  an  ideal  laser,  is  transmitted  for  each  ‘1’,  and  nothing  is  transmitted  for  each  ‘O’. 
Symbols  occur  at  regular  intervals,  ‘l’s  with  prior  probabi,:ty  f j ; ‘0’s  with  prior  probability 
ft) + £ 1 = 1*  1"  ^e  absence  of  background  noise,  we  found  that  the  signal  induces  a pure 
coherent  state  Ip)  in  the  matched  mode  of  the  aperture  field  of  the  receiver;  with  p being 
the  complex  amplitude  of  the  signal.  When  thermal  noise  of  effective  absolute  temperature 
T is  also  present,  the  density  operator  of  the  matched  mode  under  hypothesis  Hj  (signal 
present)  is 

P 1 = POO  - (1  - v)  exp  l-w(a+  - p*  1)  (a  - . (37) 

Here,  a and  a+  are  the  photon  annihilation  and  creation  operators,  I is  the  identity  operator 
in  the  Hilbert  space  for  the  quantum  states  of  the  modal  field,  and* 

w = hi>/KT,  v = e-w, 


- 


..  . . . 


i 


11 


h = Planck’s  constant,  v = signal  carrier  frequency,  and  K = Boltzmann’s  constant  The  mean 
number  of  photons  supplied  by  the  signal  is  S = Ipl  2,  and  the  mean  number  supplied  by 
the  background  noise  is  N = v/(l-v).  Without  loss  of  generality,  we  take  the  amplitude  p 
as  real  and  positive.  Under  hypothesis  Hq  (signal  absent),  the  density  operator  of  the  mode 
is 

Pg  = p(0)  = ( 1 - v)  exp  (-wa+a),  v = e-w.  (38) 

The  effective  SNR  D2  is  defined  as 

D2  = 4S/(2N+1);  (39) 

which,  in  the  limit  w < 1,  v = 1,  becomes  the  usual  SNR  2E/KT  in  terms  of  the  energy  E of 
the  received  signal. 

The  optimum  detection  strategy  and  the  minimum  probability  Pe  of  error  again, 
depends  on  the  eigenvalues  tjk  and  the  eigenvectors  1 17^  of  the  operator  p j - Xpg,  X = fg/f  1 > 

(p  J - Xpg)  I Tjk>  = T?k  1 17k>,  (40) 

as  was  shown  in  the  previous  subsections.  The  optimum  quantum  receiver  measures  the 
projector 

"1=1  lt?kXqkl,  (41) 

r,k>0 

choosing  hypothesis  Hm,  m = 0,  1,  when  the  outcome  of  the  measurement  is  m (ref  1).  The 
minimum  probability  of  error  is 

Pe  = ?lO-S  %)  <42> 

»7k>0 

or,  P^fod+X"1  I »?k).  (43) 

Tjk<0 

It  is  convenient  to  label  the  negative  eigenvalues  with  negative  indices.  When  fg  = 
fj  = Vi,  they  occur  in  pairs  of  opposite  sign  and  we  can  put  r?k  = -rj_k.  This  problem  can 
be  solved  exactly,  in  both  the  extreme  quantum  limit  (v  = N = 0)  and  the  classical  limit 
(v=  1).  For  = f,  = Yi, 

Pe  = Vi  ( 1 - ( 1 - e~S)^}  , v = 0,  (44) 

and 

OO 

Pe  = erfc  (D/2)  = (2^/  exp  (-t2/2)  dt,  v * 1 (45) 

D/2 

We  plot  our  results  in  terms  of  the  parameter 


y(v,  D)  = D2  j-41n[4Pe(l-Pe)]j  ~l, 


1 

which  is  the  ratio  of  the  SNR  D2  required  for  a given  probability  P of  error,  to  that  required 
when  no  noise  is  present  (v  = 0).  Thus,  y(0,  D)  = 1 , and 

y(l,  0)  = ir/2  <y(l,  D)  <y(l,  oo)  = 2. 

Figure  1 displays  the  parameter  y(v,  D)  in  decibels  as  a function  of  D for  various 
values  of  v,  with  dashed  lines  indicating  regions  where  our  results  lack  precision.  Figure  2 
displays  y(v,  0)  and  y(v,  «),  again  in  decibels,  as  functions  of  v = N/(N+1 ).  In  terms  of  D 
and  y(v,  D),  the  probability  of  error  is,  for  £q  = f j = Yi, 

Pe  = Vi  I1  - [ 1 - exp  (-D2/4y)J  ^ } ; (47) 

and  for  large  values  of  the  SNR  D2,  it  is  asymptotically  proportional  to  exp  [-%D2/y(v,  oo)j . 

The  error  probability  so  calculated  applied  not  only  to  on-off  signaling,  but  also  to 
decisions  between  signals  of  the  same  form  and  arbitrary  complex  amplitudes  pQ  and  p,, 
with  p = I p | - Hq\  . The  new  density  operators  are  related  to  those  in  equations  (37)  and 
(38)  through  the  unitary  transformation, 

P l'  = D0iO)p j D+(jiq),  pq  = D(p0)p0D+(p0);  (48) 


where,  D(p0)  is  a displacement  operator  whose  form  is  given  in  section  3 of  reference  1. 
Antipodal  coherent  signals  are  included  in  this  class,  with  pQ  = -p/2.  Furthermore,  when 
antipodal  two-photon-coherent-state  (TCS)  signals,  as  defined  by  Yuen  and  Shapiro  (ref  12), 
are  received  in  the  absence  of  background  noise,  the  density  operators  under  the  two  hypo- 
theses are  related  to  those  in  equation  (48),  with  pq  = -p/2  by  a second  unitary  transforma- 
tion. A unitary  transformation  of  both  density  operators  does  not  change  the  spectrum  of 
eigenvalues  in  equation  (41)  and,  hence,  leaves  the  probability  Pe  of  error  the  same. 

In  the  limit  p -*•  0 the  density  operator  under  hypothesis  Hj  is  approximately 

p 1 = p0  + Vi(poL  + Lp0)p,  (49) 

where 

L = 2(a+  + a)/(2N  + 1 ) (50) 

is  the  symmetrized  logarithmic  derivative  of  p(p)  with  respect  to  p,  evaluated  at  p = 0.  In 
the  number  representation  based  on  the  eigenvectors  I n>  of  the  number  operator  a+a,  the 
operator  p j - Xpg  whose  eigenvalues  we  seek  becomes  the  approximate  matrix, 

<n  I (p  j - Xpq) I m)  s ( 1 - X)en6nm 

+ p(2N+  ir*<nll(a+-*-a)po  + Po(a+  + a)]  lm>  (51) 

with 

e„  = (l-v)vn  (S2) 


12.  Yuen  HP,  and  Shapiro  JH,  “Optical  Communication  with  Two-Photon  Coherent  States  — Part  I: 
Quantum  State  Propagation  and  Quantum  Noise  Reduction,”  Working  Papers. 
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Figure  2.  Detection  probability  parameter  y(v,  D)  in  dB  for  v = 0 and  v ■ 1 


the  nth  eigenvalue  of  pQ.  This  matrix  is  tridiagonal;  its  diagonal  elements  are  (1  - X)en  and 
the  elements  of  the  adjacent  subdiagonals  are  (ref  2), 


<nl(p|  -Xp0)ln  + l>ip(l  - v)2  (n  + 1)*V\  n>0, 
<n I (p j -Xppjln-  l>  = /i(l  -v^n^V1"1,  n<0. 


(53) 

(54) 


(55) 


The  eigenvalues  of  a tridiagonal  matrix  can  be  calculated  by  iterative  evaluation  of  a con- 
tinued fraction  into  which  the  determinantal  equation  can  be  converted  (ref  13). 

In  particular,  when  = = = 1,  the  eigenvalues  are 

r)k  = M(l  -v)2$k; 

with  (k  the  eigenvalues  of  the  tridiagonal  matrix  I I My  I I whose  elements  are 

Mn,n+l-Mn+l.n*<"»»V, 

My  = 0,  li  -jl  # 1. 

From  equation  (43), 


(56) 
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Tnm<t»)  = WlPo  ~ Po%iO>  j ~ ’tiT^O*4!  1 m> 

= en«nm-W’J).  (60) 

with 

Wn)  = (enem)^ 1 p 1 (P 1 - »?  U~ 1 1 m> 

oo 

* (enem^  £ e^e,.- »?)■*<  nlwrXwrlm>.  (61) 

r=0 

Here,  en  is  given  by  equation  (52)  and  the  I w,)  are  the  eigenvectors  of  the  density  operator 
p,;  that  is, 

p,lwI>  = erlwr> 


er  = ( 1 - v)vr,  I wr  > = D((/)  I r >,  (62) 

in  which 

D(fi)  = exp  (/ia+  - ft*  a) 

is  the  displacement  operator.  Equation  (62)  follows  from 

p,  = D(#i)p0D+(p).  (63) 

Furthermore,  with  S = l/il  2, 

< n I wr  > = < n I D(p)  I r > 

= (n!/r!)Vi(-M*)r-ne~S/2Ln  (r‘n>(S) 

= (r!/n!)*V~re"S/2Ln  (n“r>(S) . 

In  terms  of  the  associated  Laguerre  polynomials,  a result  can  be  derived  by  expanding  the 
coherent-state  representation  of  the  displacement  operator 

D(o'^)  * exp  V^\a\ 2 + l0l 2)  <al  D(p)l/3  > 

= exp  (a*0  + or  V ~ P*0  - V4 1 p 1 2) . (64) 

In  a double  power  series  in  a*  and  0,  and  using  the  formula  of  Mollow  and  Glauber  (ref  14), 

OO  OO 

D(«*^)*  l l (n!r!)-Vn0r<nlD(p)lr>.  (65) 

n=0  r=0 

14.  Mollow  B,  and  Glauber  RJ,  Quantum  Theory  of  Parametric  Amplification,  Phys  Rev,  vol  160,  p 1076- 
1096,  25  August  1967. 


The  matrix  elements  can  be  calculated  recursively  by  means  of  the  recurrent  relations  for 
the  associated  Laguerre  polynomials. 

When  the  signal  is  strong  S > 1 the  second  term  in  equation  (60),  proportional  to 
e-^,  is  much  smaller  than  the  first;  and,  we  can  expect  the  negative  eigenvalues  of  equations 
(41)  and  (5 1)  to  be  approximately 

T?_k  as  - Xek  = - X(  1 - v)vk.  (66) 


This  approximation  can  be  steadily  improved  in  the  following  way.  For  the  eigen- 
value of  interest,  say  rj_jc,  we  substitute  - Xe^  into  the  matrix  T(rj)  of  equation  (60),  trun- 
cated to  a manageable  but  sufficiently  large  number  of  rows  and  columns;  then,  we  use  a 
standard  computer  algorithm  to  calculate  the  eigenvalues  of  the  matrix  T(-  Xe^).  The  out- 
put from  this  computation  will  be  a list  of  eigenvalues'?}^),  j = 1 , 2, . . ~ M,  where  M is  the 
dimension  of  the  matrix  and 


!>...>  |>0. 


The  kth  of  these  eigenvalues'?^),  when  thus  arranged  by  descending  magnitude,  is  a second 
approximation  to  It  is  used  to  generate  a revised  list  of  eigenvalues  by  diagonalizing 
the  matrix  Tfr/"^);  and  from  this  list,  the  kth  in  order  of  descending  magnitude  is  again 
selected.  This  procedure  is  reiterated,  until  the  value  of  the  approximation  to  ceases  to 
change  significantly.  The  eigenvalues  r?_j , j - 1 , 2, ....  M,  resulting  from  this  procedure  are 
substituted  into  equation  (43)  to  determine  the  probability  Pe  of  error. 

Diagonalizing  the  matrix  < n I (p  j — Xpg)  I m ),  representing  the  operator  in  equation 
(41)  yields  as  many  eigenvalues  as  the  number  M of  rows  and  columns  to  which  it  is  trun- 
cated. Some  of  these  will  be  positive  - half,  when  X = 1 , and  the  rest  negative.  However, 
the  iterative  method  just  described,  by  determining  M negative  characteristic  roots,  requires 
less  computing  time  for  evaluating  the  probability  of  error  to  a given  accuracy.  We  have 
found  that  this  technique  converges  quickly,  even  for  weak  input  signals.  In  addition,  it 
does  not  require  the  extensive  algebraic  manipulations  involved  in  Yoshitani’s  and  Klemm’s 
approaches  (refs  2-4).  The  solid  curves  in  the  left-hand  side  of  figure  1 were  obtained  in 
this  way. 


Because  the  matrix  elements  Snm(rj)  of  the  operator  Pg^Pj(Pi  - Tri)"*pg^  •PPe*r* 
ing  on  the  left-hand  side  of  equation  (59)  are  proportional  to  e“®,  they  will  be  much  smaller 
than  the  eigenvalues  en  of  the  density  operator  pg  when  S > 1.  The  former  operator  can 
therefore  be  treated  as  a perturbation  of  pg  in  equation  (59);  and,  as  the  eigenvalues  and 
eigenvectors  of  pg  are  known, 


Pgln  ) = enln),  en  = (1  - v)vn,  n > °,  £ en  = *»  (*7) 

n*0 

standard  quantum-mechanical  perturbation  theory  can  be  applied  (ref  15).  It  yieldt  for  the 
eigenvalues,  through  second  order  terms,  the  approximation 


15.  Dirac  PAM,  Quantum  Mechanics,  3rd  ed,  p 168-172,  Clarendon  Press,  1947. 
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^-n  en  ~ ^nn  + ^ 


'W'J-n)12 


'-ii  n nn  - — n'  ~ e - e 

m=£n  cn  em 

with  Snm  W ®ven  by  equation  (6 1 ).  The  average  probability  Pe  of  error  is  now,  by 
equation  (43),  approximately 

OO 

pe  = f0  I (« n + ^H-it) 
n=0 

r*>  * c /_  * Y *Snm^-n^2 

f0  ^nn  ^-n^  4.  e - e * 

n=0  m#n  en  em  to“) 

- 

To  terms  of  first  order  we  can  use  the  first  term  of  this  expression,  replacing  tj__  with  its 
approximate  value  -A~ *en  as  in  equation  (66), 


Pe^Pe(,)(A) 


= fo  £ X 

n=0  r=0 


ener  I <n  I D(p)  I r>  I 2 


e„  + X. 


-wwfoi  i 'MM'’ * 

n=0  r=n  r * Aennen  + xef) 

n=0  r=0  (vr-n+X)  (1  +Xvr“n)r! 

X Sr_n  jLn<r"n)  (S)J  2 

| i2~8p°>vP,iw,i> 

p=0  (vP+X)(l+XvP) 

= ViXexp  [-(2N  + 1)S1 

x f o-tpo^d^)  mA 

pK)  (vP+X)(l+XvP)  P'  ,-v' 


Here,  we  have  used  the  symmetry  of  Kn  I D(/u)  I r>  1 2 in  rearranging  the  terms  of  the  sum- 
mation, and  we  have  applied  the  formula  (ref  16): 


00 

1 


n! 

(n+p)! 


L„(p)(x)  L^y)  vn 


-a-vr1 


exp 


*x+y) 

1 -V 


(xyv)-P«  ,p(a^). 


(71) 


By  summing  equation  (69)  on  a digital  computer  and  equation  (70)  on  a programmable 
calculator,  we  have  found  that  the  second-order  perturbation  term  in  equation  (69)  becomes 
much  smaller  than  the  first  for  D 8;  and,  that  the  results  of  equations  (69)  and  (70) 
agree  within  the  accuracy  of  plotting  our  graph  for  D > 10.  The  solid  curves  in  the  right- 
hand  half  of  figure  1 were  calculated  by  means  of  equations  (69)  and  (70). 


For  D ► 1 , the  argument  for  the  modified  Bessel  functions  in  equation  (70)  is  so 
significant  that  we  can  replace  these  functions  with  their  asymptotic  form. 


IpOO  ~(2*rx)”^ex, 


in  all  significant  terms  of  the  series.  The  error  probability  is  then  asymptotically  propor- 
tional to 


1 

- 

- 

exp 

-(2N+1)S  + ^ 

= exp 

D2  (1  +v) 

4 (1+v^)2 

so  that  for  D > 1 , our  parameter  y(v,  D)  in  equation  (46)  becomes 

y(v,  “)  = (1  + v^)2(l  +V)"1  . 

This  is  plotted  in  dB  versus  v in  figure  2. 

For  large  SNR  D one  could  determine  the  false-alarm  probability  Qn  and  the  detec- 
tion probability  Qd  as  functions  of  the  parameter  X,  by  using  the  approximation  pJ^(X) 
in  the  formulas  e 


dPe  dP. 

Q0  = Pe+(l  + X)  dX’Qd=  l-pe+X(l  +x>  dx- 

These  could  also  be  applied  in  conjunction  with  previously  described  methods,  but  the 
amount  of  computation  required  would,  except  for  small  values  of  the  parameter  v,  be 
prohibitive. 


16.  Helstrom  CW,  Optimum  Quantum  Decision  between  a Pure  and  Mixed  State,  with  Application  of  TCS 
Signals,  Working  Papers. 
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DETECTION  OF  ANTIPODAL  TCS  SIGNALS 


In  order  to  transmit  a TCS  signal,  the  radiating  mode  of  the  field  in  the  transmitting 
aperture  is  placed  in  a right-eigenstate  of  the  non-self-adjoint  operator  p0a  + »>ga+;  where, 
a and  a+  are  the  usual  photon  annihilation  and  creation  operators,  and  pg  and  vq  are 
parameters  related  by  (ref  12) 

•m0  |2  - '"o  ,2  = '• 

Here,  we  assume  that  pg  and  vq  are  real.  When  such  a signal  arrives  at  a distant  receiver, 
in  the  absence  of  background  light,  it  places  the  matched  mode  of  the  receiver  aperture 
into  a mixed  quantum  state  whose  density  operator  is  a particular  unitary  transformation 
of  the  density  operator  p|  in  equation  (37)  (ref  16).  The  parameters  p and  v = e-w  of  this 
density  operator  are  related  to  those  of  the  TCS  signal  by  the  following  formulas,  in  which 
N-p  is  the  average  total  number  of  photons  in  the  received  signal  and  f is  the  fraction  of 
energy  in  the  core  of  the  transmitted  TCS  signal: 

v * N/(N+1), 


N = Vi)[4x'2NT(l-0+  ll‘A-  l}, 
p = (NT0^(m+n), 


m 


NT(l-0  + N+  1 


2N  + 1 


Vi 


, n = (m2-  1)^. 


- •* 

Here,  k1  is  the  fraction  of  the  signal  energy  lost  during  transmission  to  the  receiver,  whose 
aperture  intercepts  only  a part  of  the  transmitted  field.  In  reference  16  our  p,  m,  n,  and  N 
were  denoted  by  6,  p,  v,  and  N respectively.  Setting  f = 1 corresponds  to  transmitting  a pure 
coherent  signal  by  putting  the  radiating  mode  of  the  transmitter  aperture  in  a pure  coherent 
state  lot),  i.e.,  transmitting  a right-eigenstate  of  the  annihilation  operator  a.  The  matched 
mode  of  the  receiver  aperture  is  then,  in  the  absence  of  background  light,  also  in  a pure 
coherent  state  lp>,  and  N - 0.  For  0 < f < 1 , what  we  called  the  nimbus  of  the  transmitted 
TCS  signal  partially  deteriorates  during  transmission  into  noise,  whose  strength  is  represented 
by  the  parameter  N.  Close  to  the  transmitter,  where  this  deterioration  is  not  excessive,  the 
probability  Pe  of  error  in  detecting  a TCS  signal  can  be  significantly  smaller  than  that  of  an 
ordinary  coherent  signal  conveying  the  same  number  N-p  of  photons,  provided  the  fraction 
f is  properly  chosen  (refs  1 2,  1 5). 


Antipodal  TCS  signals  induce  two  density  operators  at  the  receiver  aperture  that 
are  the  same  unitary  transformation  of  pj  in  equation  (37),  but  with  complex  amplitudes 
+p  and  -p,  respectively.  We  assume  that  the  signals  are  transmitted  with  equal  prior  proba- 
bilities fg  ® f j = V4.  The  average  probability  of  error  will  then  be  the  same  as  what  was 
calculated  in  the  foregoing  sections,  except  that  we  replace  S * Ip  I2  by  4lp  I2  for  the  mean 
number  of  photons  in  the  equivalent  received  coherent  signal. 

In  order  to  determine  how  this  error  probability  Pe  depends  on  the  fraction  f of 
energy  in  the  core  of  the  TCS  signal,  we  have  calculated  Pe  from  equation  (70)  for  various 
values  of  the  fraction  k'2  of  energy  lost  during  transmission.  The  mean  number  N-p  of 
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photons  received  was  set  equal  to  S,  and  the  resulting  values  of  the  equivalent  SNR  D were, 
in  all  cases,  large  enough  so  that  equation  (70)  is  a good  approximation.  The  values  of  the 
error  probability  for  these  antipodal  TCS  signals  are  plotted  in  figure  3.  It  is  apparent  from 
the  figure  that  transmitting  antipodal  TCS  signals,  rather  than  ordinary  coherent  signals 
(f  = 1),  is  advantageous  only  when  less  than  about  six-tenths  of  the  signal  energy  is  lost 
during  transmission.  The  same  conclusion  w reached  in  the  limit  of  large  SNR  for  which 
equation  (70)  is  valid. 


f 

Figure  3.  Probability  Pe  of  error  in  detecting  antipodal  TCS  signals  of  prior 
probability  in  absence  of  background,  versus  the  fraction  f of  energy  in  the 
core  of  the  signals;  Ny  = 5.  Curves  are  labeled  with  the  fraction  a’^of  energy 
lost  in  transmission. 


CONCLUSION 


P 


Methods  were  presented  for  calculating  the  minimum  attainable  probability  of  error 
in  detecting  a coherent  quantum  signal  received  amid  thermal  noise  in  the  limits  of  low  and 
high  signal-to-noise  ratio.  In  the  latter  limit,  quantum-mechanical  perturbation  theory  was 
applied  to  solving  the  detection  operator  equation  approximately.  Graphic  results  are  fur- 
nished for  binary  signals  transmitted  with  equal  prior  probabilities.  Thus,  previous  work  on 
this  problem  is  extended  to  a broader  range  of  signal  and  noise  strengths.  The  results  were 
then  applied  to  the  detection  of  antipodal  two-photon-coherent-state  signals,  and  compared 
with  the  previous  results.  This  comparison  indicates  that  transmitting  antipodal  TCS  signals, 
rather  than  ordinary  coherent  signals,  are  advantageous  only  if  less  than  six-tenths  of  the 
signal  energy  is  lost  during  transmission. 
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